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Constructing high-fidelity control fields that are robust to control, system, and/or 
surrounding environment uncertainties is a crucial objective for quantum information 
processing. Using the two-state Landau-Zener model for illustrative simulations of a 
controlled qubit, we generate optimal controls for tt/2- and vr-pulses, and investigate 
their inherent robustness to uncertainty in the magnitude of the drift Hamiltonian. 
Next, we construct a quantum-control protocol to improve system-drift robustness 
by combining environment-decoupling pulse criteria and optimal control theory for 
unitary operations. By perturbatively expanding the unitary time-evolution opera- 
tor for an open quantum system, previous analysis of environment-decoupling control 
pulses has calculated explicit control-field criteria to suppress environment-induced 
errors up to (but not including) third order from 7r/2- and vr-pulses. We systemati- 
cally integrate this criteria with optimal control theory, incorporating an estimate of 
the uncertain parameter, to produce improvements in gate fidelity and robustness, 
demonstrated via a numerical example based on double quantum dot qubits. For the 
qubit model used in this work, post facto analysis of the resulting controls suggests 
that realistic control-field fluctuations and noise may contribute just as significantly 
to gate errors as system and environment fluctuations. 
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I. INTRODUCTION 



Demanding requirements for gate fidelities to achieve fault-tolerant quantum computation 
(QC) [1] have motivated the need for improved quantum control protocols (QCPs). In 
quantum information science [2], there are (at least) three distinct dynamical approaches 
to improving the fidelity of qubit operations in the presence of environmental interactions: 
dynamical-decoupling (DD) pulse sequences [3, 4], optimal control theory (OCT) [5, 6], and 
quantum error-correcting codes (QECCs) [1]. Although interesting non-dynamical methods 
exist for noise suppression, e.g., decoherence-free subspaces (DFSs) [ ], noiseless subsystems 
[8], and anyonic/topological systems [*^], the work reported in this article focuses exclusively 
on dynamical approaches for controlling quantum systems [10]. Specifically, the objective of 
this work is to construct a hybrid QCP, combining methods and results from DD and OCT, 
to locate control fields that (a) produce high-fidelity rotations about a particular axis and 
(b) are robust with respect to an uncertain frequency of rotation about an orthogonal axis. 
By combining complementary features of these analytically- and numerically-based QCPs, 
we have developed a hybrid QCP where estimates of system parameters can be directly 
incorporated into numerical simulations to generate improved quantum operations. 

During the past few years, P. Karbach, S. Pasini, G. Uhrig, and colleagues have made 
significant contributions toward the mathematical analysis and design of DD pulses and 
sequences for controlling qubit systems and decoupling them from their surrounding envi- 
ronment, e.g., refs. [1, 11-17]. In a recent article [i .], using a rather general open-system, 
time-dependent Hamiltonian for one qubit, they derive analytical control-field criteria for 
tt/2- and vr-pulses, which, when satisfied, eliminate the first- and second-order errors in 
the unitary time-evolution operator resulting from qubit-environment interactions. In this 
work, we refer to these criteria as "decoupling-pulse criteria" (DPC) and to the control 
fields that satisfy this criteria as "decoupling pulses" (DPs). We adapt the DPC for the 
case of closed-system unitary control, where the dynamics are influenced by an uncertain 
drift (i.e., time-independent) term in the qubit Hamiltonian. For control fields that satisfy 
the DPC, this adaptation eliminates the first- and second-order effects resulting from the 
drift term. Using a novel method for multi-objective control, we combine the mathematical 
DPC with a numerical procedure based on OCT for unitary control [1 8] that incorporates 
an estimate of the drift-term magnitude (i.e., system information) to construct control fields 
with increased fidelity and robustness to uncertainty in the drift term. For brevity, we refer 
to this combination of the DPC and OCT as "DPC+OCT". To demonstrate the utility of 
our approach, we optimize and evaluate these control fields using a qubit model based on 
the two-level Landau-Zener Hamiltonian [10] that has an uncertain drift term and is driven 
by a deterministic control field. Even though the qubit model is quite general, i.e., the 
Hamiltonian employed represents a one-qubit system with a linear drift term, driven by a 
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scalar control field, and thus describes a variety of qubits (e.g., atomic, spin, superconduct- 
ing, etc. [ ■]), we select physical units for the model that are relevant to double quantum 
dot (DQD) qubits to investigate the practical features of our results [20]. With this model 
system, we demonstrate that the DPC+OCT combination can be used to produce (a) im- 
proved fidelity compared to DPs alone and (b) improved robustness to uncertainty in the 
drift magnitude compared to results from DPs and OCT alone. Although research examin- 
ing the effects of classical control noise is extremely important for practical QC [21, 22], all 
control fields in this work are assumed to be deterministic, with control amplitudes that are 
exact to numerical precision. Thus, uncertainty is assumed to be present only in the system, 
and, unless otherwise specified, control robustness refers to robustness to drift uncertainty. 

Related research on hybrid QCPs includes that performed by Lidar et ai, who proposed 
the application of DD pulse sequences on logical qubits encoded in DFSs or QECCs to 
eliminate decoherence in solid-state and trapped-ion qubits [23-26]. Borneman et al. used 
OCT to design control fields that are robust to systematic amplitude and resonance inho- 
mogeneities, thereby improving the performance of the so-called Carr-Purcell-Meiboom-Gill 
pulse sequence [27]. In addition, there have been many studies on the control and control- 
lability of (inhomogeneous) quantum-mechanical ensembles, such as a collection of coupled 
or uncoupled spin systems, primarily for state-based objectives, e.g., refs. [28-32], and se- 
quences of unitary time-evolution operators that compensate for systematic off-resonant 
effects, e.g., [33-37]. 

This article is organized as follows. Section 11 introduces and develops the model qubit 
system based on the Landau-Zener Hamiltonian [ !')] used in our optimizations and simu- 
lations. For illustrative purposes, this model is compared to a logical DQD semiconductor 
qubit, where uncertainty in the drift term of the system Hamiltonian is due to the sur- 
rounding nuclear spin environment [38, 39]. A system of scaled units is defined that allows 
for the comparison of our model and control fields to relevant experimental parameters. In 
section III, we summarize our gradient-based OCT routine for deterministic Hamiltonian 
systems, describing our objective functional for unitary control, relevant control properties, 
and the numerical optimization procedure. Section IV presents results from unitary OCT 
for subsequent comparison to those from the DPC and DPC+OCT QCPs. Inherent robust- 
ness of these optimal controls (OCs) to variations in the magnitude of the uncertain drift 
term is also analyzed, and a functional is proposed to quantify this robustness. For the 
individual unitary targets considered, despite the similar structures of the resulting OCs for 
different drift-term magnitudes, their gate distances as a function of the drift magnitude 
differ dramatically. Section V summarizes the nonlinear control-field criteria developed by 
Pasini et al. for designing control pulses that are robust to decoherence [15]. Our adaptation 
of these criteria to closed-system unitary control is explained and the hybrid DPC+OCT 
control problem is posed. Control fields satisfying this criteria are applied to our model 
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qubit system and their robustness is analyzed. In section VI, we describe and mathemat- 
ically formulate our gradient-based method for solving the nonlinearly-constrained control 
problem. Section VII presents results from our DPC+OCT optimization algorithm. Gate 
distance and robustness of the control fields are numerically analyzed and discussed. We 
also compare OCT, DPC, and DPC+OCT results collectively to illustrate the benefit of 
our hybrid approach. In addition to comparing the gate distances directly, we apply these 
controls to an inhomogeneous ensemble of systems to emphasize the improvement that may 
be obtained from this hybrid QCP. Like the OCT results, significant gate-distance sensitiv- 
ities to relatively small control-field differences are observed for the DPC+OCT controls, 
supporting further study and suppression of the effects of undesired control-field fluctuations 
and noise on quantum information processing. We conclude this article in section VIII with 
a summary of our results and identify several future directions of our research. 



II. LANDAU-ZENER MODEL SYSTEM 

A. Model Hamiltonian 

We represent the dynamical model of a qubit with the following Hamiltonian (where 
h= 1; details regarding units appear in section II C): 

Hit) ■.= eS, + C{t)S,, (1) 

where Sx := (Ta/2 is a spin operator for a spin-1/2 particle, ax is a Pauli matrix (A G {x, y, z}), 
eSx represents a persistent rotation about the x-axis, and C{t) represents the time- dependent 
control field driving rotations about the 2;-axis. Note that H(t) corresponds to the two-state 
Landau-Zener model [19], and that both C(t)/h and e/h have units of angular frequency, 
e.g., radians per second in SI units. 

Let T-L denote the Hilbert space of the system, where n := dim{H} [n = 2 for one qubit), 
and denote the orthonormal basis of Sz that spans "H, with corresponding eigenvalues 

±1/2. The Lie group of all unitary operators on "H is denoted by U('H). In general, the 
unitary time-evolution operator U{t) G U('H) for a closed quantum system obeys the time- 
dependent Schrodinger equation: 

U{t) = -iH{t)U{t), (2) 

where U{t = 0) = ll„, the n x n identity matrix. From a controllability perspective [40, 41], 
the Hamiltonian in eq. (1) generates the Lie algebra su(2). Thus, the system is completely 
dynamically controllable, i.e., any element of the Lie group SU(2) can be generated via 
eq. (2) and an appropriately-shaped control field. However, this analysis does not necessar- 
ily reveal anything about the control-field structure required to realize an arbitrary SU(2) 
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operation. As an illustrative example of our DPC+OCT QCP, we focus on constructing 
unitary operations corresponding to 7r/2- and vr-rotations about the z-axis. 

B. Double quantum dot logical qubit 

Although our qubit model is quite general, for illustrative purposes we refer to a particular 
application of a DQD solid-state qubit [42, 43], which has been studied in an array of 
experiments, e.g., [44-48]. With one electron in each quantum dot, the DQD system spans 
four spin-1/2 states. An applied magnetic field will break the degeneracy of the states in 
which both electrons are either aligned against or with the field. In this situation, it is 
possible (and often advantageous) to work within the two-level subspace where the net spin 
angular momentum is zero. By adjusting voltages in the electrostatically defined quantum 
dots, the magnitude of the exchange interaction between the electrons may be controlled. 
This interaction controls the splitting between the singlet, \S) := (|ti) ~ IJ't))/"\/2 and 
triplet, |To) := (|ti) + \it))/V2, states of the spin-zero manifold. Designating \Sq) = \S) 
and |iSi) = \Tq), we equate C(t) in our general model [eq. 1] with the exchange interaction. 

The spin-zero manifold is insensitive to a global magnetic field. However, gradients in 
the magnetic field will cause singlet-triplet transitions. Such a gradient splits the energy of 
the states \U) = (l^o) + \S))/V2 and |;t) = (l^o) - \S))/V2 by the difference in effective 
Zeeman energies for an electron in either of the two quantum dots. In this context, we may 
therefore equate the energy e with this effective Zeeman energy difference. In GaAs DQD 
systems, the effective Zeeman energy difference is typically dominated by the Overhauser 
shifts from a lattice of randomly polarized nuclear spins corresponding to approximately 
1.602 X 10-26 - 1.602 X 10-25 J (or 10-^ - 10-^ eV) [42, 43]. It has been demonstrated that 
a desired difference in Overhauser shift of a GaAs DQD may be realized through feedback 
control from a preparatory qubit [ ]. However, the value of e will drift over time through 
the nuclear spin diffusion that causes spectral diffusion [3, 38, 39, 49], motivating the need for 
robust control. In proposed Si DQD systems, e.g., ref. [oO], nuclear spins may be eliminated 
through isotopic enrichment. Other spin baths, such as electron spins of donor impurities 
[51] or dangling bond spins at an interface [52], may also lead to variations and drift in the 
value of e. 



C. Scaled unit system 

In addition to setting the reduced Planck constant h = 1 (corresponding to unit of 
angular momentum: energy x time), a simple set of scaled units is defined by also setting 
the final time of the controlled evolution tf = 1. The ratio h/t{ = 1 yields the scaled unit of 
energy, which, when tf = 20 ns (as an example), corresponds to approximately 5.273 x 10^^^ 
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J (or 3.291 X 10~^ eV). By appropriately scaling the Hamiltonian and control field, this 
propagation time can be transformed to any final time tf. Relationships between these 
scaled units and SI units (when tf = 20 ns) are summarized in table 1. 

In this work, optimizations were performed for individual values of e G [0,5]; we denote 
the nominal values of e used in these calculations as Eq- This range of e corresponds to no 
and moderate rotations from the environment for the GaAS DQD example. For a DQD 
logical qubit, 

e = gelJ'BBA, (3) 

where ge is the so-called electron (^-factor, /ib is the Bohr magneton, and is the magnetic 
field resulting from the difference in the random hyperfine fields from each quantum dot 
along the direction of the applied field. When 1 scaled unit of time corresponds to 20 ns 
(a representative estimate of the time required for one-qubit rotations for a DQD system 
[45, 53]), e = 5 scaled units of angular frequency (the maximum value of considered) 
corresponds to Ba ~ 6.5 mT; this is consistent with experimental reports of GaAs DQDs 
(where ge = —0.44) [42, 43, 45]. Unless stated otherwise, all physical quantities in this work 
are expressed in scaled units. 



Physical quantity 


Scaled unit 


SI unit 


angular momentum: h 


1 


1.055 X 10-34 Js 


time: t 


1 


2.0 X 10"^ s 


energy: C, e 


1 


5.273 X 10-27 J 



TABLE I: Scaled and SI units for the logical qubit described by the Hamiltonian in eq. (1). 



III. OPTIMAL CONTROL OF UNITARY OPERATIONS VIA 
GRADIENT-BASED ALGORITHMS 

A. Objective functionals for unitary operations 

For a target unitary operation V G U(?^), the distance A between V and a simulated 
final-time unitary operation U{t{) is 

A[V, UiU; C)] := min ^ ||f/(tf; C) - exp(i^)y||Hs , (4a) 



l--|Tr[Vtf/(tf)]|, (4b) 



n 



where || ■ ||hs denotes the norm based on the Hilbert-Schmidt inner product: (A, B)hs '■= 
Tt{A^B), A,B E A^„(C) [A^„ (X) denotes the set of n x n matrices over the field X]. 



7 



This phase-invariant distance measure is a special case of a more general distance measure 
developed in ref. [ > 1], which is applicable to studies involving composite systems where only 
the qubit/system dynamics are directly of interest [55, 56]. 

Concerning mathematical notation, because the unitary time-evolution operator is a func- 
tion of time and a functional of the control, it will be expressed more generally as U {t; C) 
for all time t and a control C, compared to U{t)\ the final-time unitary operator will be 
expressed more generally as Utf{C), compared to U(t{). Also, we denote the space of ad- 
missible controls with final time t = tf as €tf Some properties of the Hilbert space (tfj are 
discussed below; further details are in ref. [ > i]. 

Because < A < 1 in general, it is useful to define the fidelity J-" of unitary operations 
as [54, 57] 

^ := ^ |Tr [0t/(tf)] 1 = 1- A\V, t/J, (5) 

which is a common phase-invariant measure of gate fidelity based on the Hilbert-Schmidt 
inner product, e.g., refs. [18, 58, 59]. Note the quadratic dependence of J-" on A, i.e., a 
distance of 10"^' corresponds to a fidelity of 1 — 10^^, where x < 0. 

An optimal control field for a given unitary operation may be located by minimizing an 
objective functional J'[C] of the control field that incorporates the final-time unitary target 
V, constrains the dynamics of U(t) to evolve according to eq. (2), and penalizes the fiuence 
of the control field. For this work, the objective functional is defined as 

J[C] := A[V, U,,iC)] + I l^'^dt. (6) 

Often, the minimization of is performed using a gradient-based algorithm (GrA; see 
[18, 54, 55, 60] for details and examples of gradient-based optimizations). Here, a > weighs 
the control-field fiuence relative to the distance A and s : [0,tf] i-)- M is a continuous "shape 
function". When appropriately chosen, s{t) penalizes undesirably-shaped functions [5 I]. We 
use s{t) = sm^{nt/ti), where p G Q+ U {0}. For p ^ 0, this form penalizes the control-field 
slew rate around the initial and final times and favors controls where C(0) = C{t{) = 0. 

For the time-dependent Hamiltonian in eq. (1), Utj : Ct^ — > U(7{) denotes the map, defined 
implicitly through the Schrodinger equation [eq. (2)], that takes a control field C G to 
the final-time unitary evolution operator Utf G 11(7/). Note that is a Hilbert space of 
admissible controls, on which U{t; C) exists for all t G [^lU] and all C G [ ], where the 
inner product on (tjj is 

{f,9K-= ly-^J^^t' V/,^?G£,, (7) 

As such, JT" : — )■ M is the dynamical version of the distance measure A, with a relative 
cost on the control-field fiuence, determined by a. The role of the shape function s{t) in 
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eqs. (6) and (7) is to change the geometry of control space, moving undesirably-shaped 
functions away from the origin, out to infinity, where they are less likely to be the targets 
of a minimization over [54]. 



B. Control rotation angle 6(t; C) 

In addition to the objective functional and inner product on €tf, another important 
expression is the integral of the control field: 

0{t; C) := / C(r)dr. (8a) 
Jo 

The angle 6 corresponds to the rotation about the z-axis performed by the control field 
during the time interval [0,t]. Although ^ is a functional of the control field C(t), when- 
ever appropriate we abbreviate 9{t; C) as 9{t) to avoid unnecessarily cumbersome notation. 
Equation (8a) is equivalent to 

d9{t;C) 
dt 

Also, 

5e[t;C] 
6C{t') 

(i.e., the functional derivative of 9 with respect to C) where S^{t — t') is the Heaviside step 
function: 

, , , ,1 when t>t' 
^ ^ ^ when t<t' ^ ^ 

C. Optimization with a gradient-based algorithms 

This section briefly summarizes the variational analysis of and describes criteria for 
the optimal points (or submanifolds) of with respect to a control field C. The gradient of 
the objective functional is explicitly derived in ref. [ > l]; we present it here for continuity: 



C{t). (8b) 
^{t-t'), (8c) 



{VJ[C]){t) = 4^^[^|^^ ^^^j Im(Tr{[[/t.(C)i?-i?t^^^(C)]yt(^.C>)g^^(^.C>)|j 

(10) 

where R := exp{iLp)V and (p := Imjln [Tr (V^"''[/fj)] }. Critical points of J' (a real-valued 
functional) are defined as controls for which (VJ'[C])(t) = for all time t [(>2]. Control fields 
are iteratively updated using this gradient to improve the value of the objective functional 
J. Given the kth iterate of the control field C^''\t), adjustments to the control field for the 
(/c + l)th iteration are given by 

C'-''+^\t) := C^''\t) - 13 {VJ [C(^)]) (t), (11) 
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where /9 is a constant that determines the magnitude of the field adjustment. This procedure 
describes an implementation of a steepest-descent algorithm [63] . In this work, initial control 



final times and slew rates are consistent with the shape function s{t) = sin(7rt/tf). 

IV. RESULTS FROM QUANTUM OPTIMAL CONTROL THEORY 

Using only the GrA presented in section III, OCs were found for unitary targets that 
perform 7r/2- and vr-rotations about the z-axis: 



where G {vr/2, vr}. The final time for all OCs was fixed at tf = 1 scaled unit of time. With 
the GrA and the objective functional J', a combination of the value of e and the structure 
of the initial control field determines the resulting optimal control field. Optimizations were 
performed individually for specific angular frequencies: Eq G [0,5]. As described in section 
II C, this interval represents the regime of zero to moderate rotation from the environment 
for the DQD logical qubit. To emphasize the e-specific nature of these OCs, we denote them 
as Co{eo;t). 

Because of the similarity of results over the entire interval < e < 5, only a subset 
will be presented. OC fields for and Z^^ as a function of e are presented in figs. 1 
and 2, respectively, for Eq G {0,1,2,3,4,5}. Even though all of these OCs were located 
using the same initial control field, which is very similar to the OC reported for Eq = for 
both targets, some of the converged fields differ dramatically for different values of Eq. All 
OCs have distances A < 10~^ (or J-" > 1 — 10~^^, essentially corresponding to the limits 
of numerical precision), which is expected because this system is relatively simple and fully 
controllable [40, 41]. For Eq = 0, all OCs for satisfy 9(tf) = (j) (mod 27r), i.e., OC design 
simply corresponds to pulse-area control in this situation. However, when e ^ 0, there is 
no corresponding pulse- area requirement. In fact, if e is known accurately, it is possible 
to perform Z^ operations with piecewise constant controls that satisfy 0{t{) = 0. Table II 
contains information about some of the properties of these OCs. For a DQD logical qubit, 
we observe that these controls require negative exchange coupling values. Although negative 
exchange energy is uncommon, it is predicted to be possible to produce through combined 
tuning of the magnetic field, dot size, and tunnel coupling [64]. 

For both and operations, despite the similar structures of the OCs, especially 
Co(l; t) and Co(2; t), their gate-distance responses for < £ < 6 (with a numerical resolution 
of 0.01 scaled units) are quite unique, as shown in figs. 3 and 4, respectively. Even though 
max^ |Co(l; t) — Co{2] t)\ < 0.2 scaled units of energy for both operations (see inset of figs. 1 



fields C-^^ are continuous approximations to simple square-wave pulses, where initial and 




(12) 
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and 2), and the mean relative difference is approximately 1.4% and 5.7% for the Z^^ji and 
operations, respectively, this two-level system effectively discriminates between these two 
similar control fields through the response of the distance functional A in eq. (4). Within 
the interval < e < 3, the gate distances of Co(l; t) and Co(2; t) do not significantly overlap. 
The sensitivity of this system to these relatively small control-field variations combined with 
the inherent noise (and limited resolution) present in most realistic control sources warrants 
further study of the impact that realistic control-field fluctuations may have on practical 
fault-tolerant QC [20, 50, 58, 65]. 

For the target operation when ^ {0,1,2,3}, the OCs produce a net positive 

angle of rotation about the 2;-axis, given the initial control field. When eo ^ {4,5}, OCs 
produce a net negative angle of rotation about the 2;-axis. As the static angular frequency of 
the rotation about the x-axis increases, the OC strategy tends toward a controlled rotation 
about the z-axis in the negative direction. OC simulations for 3 < eo < 4, with a numerical 
resolution of 0.01 scaled units of energy (detailed results are not reported), reveal a distinct 
transition between OCs with shapes very similar to Co(3,t) and Co(4, t) reported in fig. 1, 
corresponding to net positive and negative rotations, respectively. For the initial control 
field used in this work, this transition occurs when Eq ~ 3.9. Comparing the gate-distance 
responses in fig. 3 for Co (3, t) and Co (4, t) reveals that these OCs are not equivalent solutions 
for unitary control; the difference between these controls is exclusively due to the effect of 
the different values of oii the qubit dynamics. 

Quantum-computing architectures often assume encoded quantum operations to correct 
the inevitable errors due to control and environmental noise. Gate operations, such as 
simulated here, must achieve a minimum fidelity threshold for successful quantum error 
correction. A predicted maximum distance A of less than 10~^ for e G [0, 5] is within typical 
ranges necessary for fault-tolerant QC [1, 43, 50]. These results highlight the potential 
advantage of using OCT if estimates of Hamiltonian parameters are known well. However, 
these OCs are not robust to uncertainty in the magnitude of e\ in fact, they are highly 
sensitive to small perturbatives in e. Uncertainty can result from incomplete or poor system 
parameter estimates as well as from dynamics of the environment [38, 51]. 

The objective functional J in eq. (6) does not include criteria to evaluate control-field 
robustness with respect to variations in e. To investigate any inherent robustness quan- 
titatively, OC fields optimized for particular values of e [i.e., Co(£o;^)] were subsequently 
applied to a surrounding interval of Eq. Results are presented in figs. 3 and 4 for the 
and operations, respectively. Consider the distance response of Co(2; t) to variations in £, 
which varies substantially with respect to variations that correspond to ~1 mT fluctuations 
(corresponding to approximately 1.6 < e < 2.4) for GaAS DQD systems. Without speci- 
fying a measure of robustness as an additional control objective, the resulting OCs are not 
inherently robust to modest local magnetic- field fluctuations (e.g., at each quantum dot). 
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Target operation: 








1 


2 


3 


4 


5 


max Co 


2.4 


18.5 


18.4 


18.3 


6.7 


5.9 


0{tf;Co) 


tt/2 


1.5779 


1.6002 


1.6406 


-3.4333 


-2.5660 




3.4 


84.5 


83.0 


80.7 


17.5 


12.1 




4.59 X 10"*^ 


1.11 X 10"'^ 


2.03 X 10-^ 


6.23 X 10"^ 


1.25 X lO-'^ 


3.94 X 10-® 


■^71/2 [Co, £0,0.5] 


7.40 X 10-2 


1.44 X 10^3 


3.65 X 10-3 


6.96 X 10-3 


5.04 X 10-2 


6.21 X 10-2 



Target operation: Z-j^ 








1 


2 


3 


4 


5 


max Co 


4.8 


12.2 


12.1 


11.9 


11.6 


11.4 


e{tf,Co) 


vr 


3.1020 


2.9824 


2.7802 


2.4923 


2.1181 




13.7 


39.5 


38.4 


36.5 


34.1 


31.4 




6.05 X 10-*^ 


2.73 X lO-'^ 


3.16 X 10-® 


5.58 X 10-® 


4.34 X 10"® 


2.79 X 10-® 


7^,[Co, £0,0.5] 


3.80 X 10-2 


6.95 X 10-3 


1.38 X 10-2 


2.05 X 10-2 


2.68 X 10-2 


3.24 X 10-2 



TABLE II: Performance of the OCs Co{eQ;t) for one-qubit operations. Here, max |Co| 
e, $[C] := J^'C\eo;t)dt, A, and are the maximum control-field amplitude, angle of 
controlled z-axis rotation, control-field fluence, gate distance, and gate robustness, 
respectively, in the corresponding scaled units described in section II C. 



Numerical calculations with these OCs for both operations indicate that the error in the 
measurement fidelity used to characterize e for a particular system must be smaller than 
10"2 (corresponding to approximately 1.3 x IQ-^ or smaller T for the GaAs DQD example) 
to realize gates with distances that are IQ-^ or smaller. Moreover, with these controls, e 
could not drift significantly during a computation without serious fidelity loss. 

In addition to the data presented in figs. 3 and 4, we introduce the following functional 
to investigate robustness of operations over the interval [£_,£_).]: 

7^^[C, 6, 6e] := ^ [Z^, Ut,{C)] de', (13) 

where e± := e ± 6e. Values for 7^0[Co, £0, 0.5], corresponding to the average gate distance 
over a unit interval centered at Eq, are reported in table II. Quantifying robustness with this 
metric further demonstrates the general lack of robustness of these OCs; 7^^ varies from 
1.44 X 10-3 (which is somewhat robust) to 7.40 x IO-2. 
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V. ROBUST DECOUPLING-PULSE CRITERIA 



Optimization of the functional J in eq. (6) is highly under-determined, and multiple con- 
trol fields exist that will produce the same target operator V [59, 66]. Requiring robustness 
to control and/or system variations involves the specification of additional constraints or 
penalties, such as eq. (13), thereby limiting solutions to this OCT problem. In this section, 
we summarize a set of control-field constraints that characterize robustness to perturbative 
decoherence and adapt them to locate controls that are robust to system uncertainty. 



A. General robustness criteria 

Consider the following open-system Hamiltonian for one qubit: 

i^open(t) := S ■ C{t) + S-f + H,, (14) 

— # — * 

where S := {Sx,Sy,Sz) represents the spin-operator vector, C(t) := (Cx,Cy,Cz) represents 
a multi-polarized control field, F := (F^., F^, F^) represents the environment interaction op- 
erator, and He represents the environment Hamiltonian. 

By expanding the final-time unitary evolution operator generated by -ffopen with respect 
to tf||iJe|| and tf||F|| about tf||ifc|| = and tf||F|| = 0, Pasini et al. have identified control- 
field criteria necessary to eliminate perturbative first- and second-order effects resulting from 
the environment Hamiltonian and the interaction term S" ■ F [ ]. Although applicable 
to multi-polarized controls and general qub it-environment coupling, the methodology devel- 
oped in this section assumes control-qubit coupling along the 2;-axis and qubit-environment 
interaction along the 

^open('^) := C{t)Sz ® tn. + ® F, + I2 ® H,, (15) 

where rze := dimj^e} and "He is the Hilbert space of the environment. For controlled 7r/2- 
and TT rotations about the 2;-axis, Pasini et al. derived the following vector functional char- 
acterizing the space of controls that suppress first- and second-order effects errors resulting 
from an x-axis interaction with the environment: 

C>]:=(Ci, (2,(3,(4,(5)^, (16a) 
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where 




■2, 



(16b) 



(16d) 



(16c) 



(16e) 



(16f) 



Recall 6{t; C) from eq. (8a) in section IIIB, which represents the net rotation performed by 
the control field C during the time interval [0,t]. 

For convenience and simplicity in the analysis that follows this section, we define r/ as 



Thus, rf : (ttj H- M^. For one qubit, the components of fj represent the first- and second-order 
perturbative errors, with respect to the final time t{ and error Hamiltonians F^; and ife, of 
a controlled 7r/2- or vr-rotation about the z-axis. Specifically, rji and 772 represent first-order 
errors, while 773, 774, and 775 represent second-order errors. Thus, when ff = 0, the pulse 
is accurate up to third-order, eliminating the first- and second-order effects resulting from 
a perturbative qubit-environment interaction. According to the analysis in ref. [1 j], when 
[Hei Fa] = 0, for all A, components 774 and 7/5 can be neglected from the vector constraint. 



To apply these results to a closed one-qubit system and construct robust operations 
for the DQD logical qubit using this criteria, we first compare the Hamiltonians H in eq. (1) 
and -ffopen ^l- (15). These Hamiltonians are equal if F^,. = e and = 0, which implies 
that [ifg; r^] = 0, so "ff^ := (771,772,7/3)"'" is the relevant reduced vector constraint. Incorpo- 
rating these nonlinear equality constraints into the original optimization problem yields the 




(17a) 



where 



r]i := (i o e. 



(17b) 



B. Closed-system robustness criteria 
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following nonlinearly-constrained control problem: 

min J\C\ 

subject to (18) 

m = 

Methods such as "diffeomorphic modulation under observable-response-preserving homo- 
topy" (DMORPH) [58, 67-69] or sequential quadratic programming [70] are required to 
generate OCs that maintain or satisfy approximate feasibility, determined by ff^ = 0. A 
technique using DMORPH is developed and applied in the next sections. 

For a qubit described by the Hamiltonian in eq. (1), control fields from ref. [ J 5] satisfying 
||57''||2 < 10^"^ (where || ■ ||2 denotes the vector two-norm) and corresponding gate-distances 
for tt/2- and vr-pulses are presented in figs. 5 and 6, respectively. These fields are denoted as 
Cd(t), where the subscript "d" indicates the decoupling feature of these DPs. Satisfying the 
control-field constraints specified by Pasini et al. [ I 5] , first- and second-order perturbations 
about e = are eliminated. As such, gate distance increases with the magnitude of e and 
optimum performance occurs when e = 0, which is not necessarily expected to be valid for 
reahstic qubit systems with drift terms, e.g., [15]. Although we use = as a general 
condition for robustness in this work, whether controls satisfying ff = are robust abouts 
points where e 7^ remains an open question. However, comparing the gate distances in 
fig. 6 to figs. 3 and 4 reveals a certain degree of robustness in control fields Cd(t) relative 
to Co(£;t), e.g., for both operations and all values of Eq considered, dA/de around Eq 
for Cd(t) is much smaller than Co(£o;^)- Table 111 contains information about some of the 
properties of the DPs in figs. 5 and 6. 



VI. HYBRID QUANTUM CONTROL: DECOUPLING-PULSE CRITERIA 

OPTIMAL CONTROL THEORY 

Given the favorable structure of quantum-control landscapes, e.g., trap- free structure, 
continua corresponding to optimal solutions, etc., for regular controls [59, 65, 66, 71], 
DMORPH provides a mathematical means to explore families of controls that achieve the 
same objective [58, 67]. Applications of DMORPH include the continuous variation of a 
Hamiltonian while preserving or optimizing the value of a quantum-mechanical observable 
[67, 72] and exploring the level sets of state and unitary control [58, 69, 73]. DMORPH 
can also be used as direct optimization technique [7 1] . We develop DMORPH techniques to 
explore the space of controls corresponding to i/ = while optimizing J for a specified Eq. 

Expressed more formally, in this section, we develop a method to optimize J over the set 
^tf •= ^ ^tf) the set of feasible controls satisfying ff = 0, where f]~^ denotes the 



15 



Target operation: 








1 


2 


3 


4 


5 


A(Z,/2,C/t,) 


7.60 X 10"*^ 


1.52 X 10"^ 


1.34 X 10"^ 


4.47 X 10-3 


1.02 X 10^2 


1.89 X 10^2 


■^7r/2[Cd, £0,0.5] 


2.84 X 10^6 


1.96 X 10"^ 


1.42 X 10"^ 


4.58 X 10-^ 


1.03 X 10~2 


1.91 X 10^2 


max |Cd| = 29.5 6l(tf ; Cd) = 7r/2 $[Cd] = 335.5 ||r?'"(t; Cd)||2 = 2.52 x 10"^ 



Target operation: Z„ 



£0 





1 


2 


3 


4 


5 




5.67 X 10-s 


5.84 X 10-^ 


4.63 X 10-3 


1.54 X 10-2 


3.56 X 10-2 


6.75 X 10-2 


7^,[Cd, £0,0.5] 


1.83 X 10"^ 


7.29 X 10-^ 


4.91 X 10-3 


1.58 X 10-2 


3.61 X 10-2 


6.81 X 10-2 




max|Cd| = 28.8 


e{tf,Cd) = 7T 


$[Cd] = 264.8 


||r/'"(i;Cd)||2 = 8.04 X 10-8 





TABLE III: Performance of the DPs Cd(t) for one-qubit operations. Here, A, TZ^, 
max |Cd|, 9, <I>[C] := f^'C\eo;t)dt, and || 7/''||2 are the gate distance, gate robustness, 
maximum control-field amplitude, angle of controlled z-axis rotation, control-field fluence, 
and constraint vector norm, respectively, in the corresponding scaled units described in 

section II C. 



preimage of ff ^. To increase general applicability, we develop this technique for ff, rather 
than the reduced constraint fj^. Away from critical points of fj (a real- valued vector), i.e., 
controls for which the set of gradients {V?7i} are linearly dependent [75], (t^^ is a codimension 
5 submanifold of (Etf [76]. It is assumed that critical points of ff are rare, an assumption 
supported by the success of the resulting algorithm. The gradient of the restricted functional 
/C := J'l^v at a point C G is just the projection of the gradient of JT" at C onto the tangent 
space Tc^ti C [58, 77, 78]. By systematically updating the control field iteratively 

using a GrA with this projected gradient, the algorithm is able to simultaneously improve 
the value of J and maintain approximate feasibility, or at least impede deviations from 
feasibility. It is unlikely that the quantum-control landscape for the restricted objective 
functional /C is trap-free. As such, a global optimization algorithm might be better suited 
to finding solutions. However, because a global parameterization of the set is lacking, 
maintaining (approximate) feasibility, i.e., ff = 0, might be difficult in general. 



The preimage of a particular subset S C Y of the codomain of a function / : X i-> F is the set of all 
elements of the domain X oi f that map to elements of 5", i.e., f^^(S) -.^ {x E X : f{x) e S*}. Because 
f] : €tf ^ M^ this implies that r/-^ -.R^ ^ 
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A. Gradients of the feasibility constraints 



Using DMORPH to remove the components of VjT that cause a change in ff requires the 
gradient of each element of ff, Vrjf. 

[''Km 
It mr 

( \ 



5C{t) 



;i9a) 



( (Vr/i[C])(t)^ 

v(v%[c])(t) J 



cos[6'(r)]dr 
sin[0(r)]dr 



2 / Y '{ cos[^(r)] cos[^(r )] + sin[0(r)] sin[^(r )] }sgn(r - r')dr'dr 
Jt Jo 



\ 



T cos[^^(r)]dr 

f 

r sin[6'(r)]dr 



(19b) 



As expressed here, the vector of gradients Vr^j is a function of the time variable t. Note that 
the set {Vrji} spans the normal space (Tc-C^J''' when C G is a regular point of fj. 



B. Gradient projection method 

In addition to the gradients Vr^j, we also need a vector that specifies the relative weight 
of each gradient component to remove. This is determined by first calculating the Gramian 
matrix, with elements 

{Gc),,--={^m[ClVr^AC])^,^- (20) 
In general, Gc is not guaranteed to be full rank; non- singularity of Gc must be explored 
(numerically) as a function of G . The Gramian matrix Gc is rank deficient if and only if 
elements in the set {Vr^j} are linearly dependent, i.e., if and only if C is a critical point of 17 
[75]. However, when Vr^j are all linearly independent, Gc is invertible. With eqs. (19) and 
(20), all gradient directions Vrji can be removed from VJ, producing V/C as follows: 

5 

V/C[C] = VJ[G] -Y^Vr^AC] {G-^ % {VJ[C])]]^ , (21a) 

i=l 

where /C : <tl^ — )■ M is a restriction of J , i.e., /C := J\tl ? and has elements 

[^c {yJ[C]% := {VJ[Gl VnAC])^,^ . (21b) 
Thus, V/C is a vector field on (tj^, and the ordinary differential equation (ODE) 

dC(s) 



ds 



-V/C[C(s)] 



(22) 
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describes the gradient flow of /C on (t^^ tliat minimizes J witliout clianging tlie value of r/. 
The GrA in this work implements a forward Euler integration of this equation, which should 
be sufficiently accurate, provided that the multipher /3 in eq. (11) is selected properly, i.e., 
/3 is within the validity of the linear approximation of the tangent space at C . Higher-order 
numerical ODE solvers, e.g., Runge-Kutta methods [63], might offer higher accuracy and/or 
greater efficiency, but have not been explored in this work. 

Equation (21) describes the orthogonal projection from TcCf^ to TcC)'^. That V/C is 
orthogonal to all elements of {Vr/i} can be verified as follows. Let 

5 

e:=5Zx.Vr74^], (23) 

i=l 

i.e., ^ is a linear combination of the elements of the set {V?7i}, where Xi ^ I^- Replacing 
SI J with ^ in eq. (21) yields 

5 5 

i - 5^ Vr/JC] {G-^ % (0]}, = i - 5^ Vr/JC] {Q-'G^x), = 0, (24) 

i=l 1=1 

where x '■= (Xi? X2, Xs, X4! X^)'^- If ^ ^ T^c^ti is such that ^ is orthogonal to all gradients Vrji, 
then the projection described in eq. (21) acts as identity on C,- Together with the previous 
statement, this shows that eq. (21) is the orthogonal projector from Tc^tf to Tc^tc 

As mentioned in the introduction, because we combine DPC and OCT methods to gen- 
erate improved control fields, we denote the integrated optimization procedure described in 
this section as DPC+OCT. Straightforward modifications of Gc and Qc are required when 
■ff^ is the constraint vector rather than ff, i.e., Gc € A^3(]R), rather than A^5(]R), and 
Qc {VJ) e ]R^ rather than 

VII. RESULTS FROM DECOUPLING-PULSE CRITERIA + OPTIMAL 

CONTROL THEORY 

Using the DPC+OCT protocol described in section VI and the DPs in fig. 5 for the 
initial iterations of all values of Eq considered, we sought to numerically explore the space 
of controls satisfying ff^ = and A = to improve control fidelity and robustness to e- 
uncertainty for and operations, compared to the original DPs. To a certain extent, 
it appears a priori that the minimization of A and ff^ might be competing control objectives. 
For example, compare the gate-distance plots for OCT (figs. 3 and 4) to those for the DPs 
(fig. 6) over the interval < £ < 6 (with a numerical resolution of 0.01 scaled units). 
OCT for the design of unitary operations, as we have formulated it in section 111, seeks to 
minimize A for a particular value of e, namely, the parameter estimate Eq. Because the 
system described by the Hamiltonian in eq. (1) is controllable and the underlying control 
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landscape possesses a fortuitous structure for regular controls [59, 66], a GrA achieves this 
objective quite efficiently and successfully. However, as presented in section IV, these OCs 
are not inherently robust to perturbations in e, whereas controls satisfying ff^ = are nearly 
optimal (with respect to J') when eo ~ 0. This scenario illustrates the potential balance 
that can exist between fidelity and robustness in general. Overall, the controls that satisfy 
= are reasonably robust to perturbations in e, provided that e and/or Ae are within 
the perturbative limit. Given that these two control objectives are potentially competing, 
we use the hybrid optimization procedure developed in section VI to suppress deviations 
from rf^ = 0, but not entirely eliminate them. In other words, convergence of the DMORPH 
DPC+OCT algorithm occurs only when J' stops decreasing, not when ff^ starts increasing. 

Because the formulation of OCT requires the specification of e, we consider Eq G [0,5]. 
As before, the final time for all controls was fixed at tf = 1 scaled unit of time. DPC+OCT 
fields for and are presented in fig. 7 and fig. 8, respectively. These fields are denoted 
as Ch(eo; t), where the subscript "h" indicates the hybrid feature of this QCP. Corresponding 
results for the gate distance A are presented in figs. 9 and 10; results for the vector-constraint 
norm ||57'^||2, and objective functional J7 are presented in figs. 11 and 12. In addition, table 
IV contains information about some of the properties of the DPC+OCT fields. For both 
operations, when Eq > 2, we note that ||57'^||2 > 10~^. However, since the robustness criteria 
quantified by ff^ were obtained from a perturbative analysis about Eq = 0, it remains an 
open question whether, a priori, control fields satisfying 17 ^ = for > will be robust 
to fiuctuations about Eq. However, as we present in this section, control solutions obtained 
from the DPC+OCT protocol have some desireable properties, even when Eo > 0. It is 
interesting to explore the results of the DPC+OCT protocol when the distance is relatively 
large (i.e., A > 10^^) and the sensitivity to changes in e is small, e.g., as illustrated in fig. 6. 

Despite their unique gate-distance dependence of these control solutions on e, as shown in 
figs. 9 and 10 over the interval < e < 6 (with a numerical resolution of 0.01 scaled units), 
the converged DPC+OCT fields for eq G {0, 1, 2, 3} are very similar to each other and the 
originated DP for each target unitary operation. This supports the observation in section IV 
that this simple system can effectively discriminate between very similar control fields, i.e., 
as measured by the gate distance A in eq. (4), the qubit system is quite sensitive to these rel- 
atively small control-field variations. For example, although maxj |Ch(l;t) — Ch(2;t)| < 0.2 
scaled units of energy for both operations, and the mean relative difference is approximately 
1.5% and 1.3% for the ^7^/2 and operations, respectively, the corresponding gate distances 
(presented in figs. 9 and 10) do not coincide significantly when < £ < 3. 

Interestingly, the DPC+OCT control field for the Z,^ operation when Eq = 5 has some 
distinguishing features compared to the fields for the other values of Eq- Based on the 
relatively large distance A of the corresponding Cd(t) control used for the initial iterate in 
the DPC+OCT protocol (A = 6.75 x 10^^), this value of Eq is not within the perturbative 
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limit of the analysis that produced the vector constraint ff = 0. With such a large distance 
at Eo = 5, the DPC+OCT routine improves the distance by a factor larger than 10^ and 
simultaneously improves the robustness TZt^ for 4.5 < e < 5.5 by a factor larger than 10, 
compared to the corresponding results for Cd(t) presented in table 111. 

Figures 9 and 10 compare the distance of Cd(t) and Ch(£:o;^) control fields for both 
operations. Compared to Cd(t), all control fields Ch(£:o; t) for > exhibit improved robust- 
ness to ^-uncertainties in an interval around the nominal value Eq used in the DPC-I-OCT 
algorithm. This result demonstrates the utility of combining so-called "pre-design" meth- 
ods, which are based on mathematically analyzing general models [e.g., eq. (14)], such as the 
DPC developed by Pasini et al. [J""], with numerical OCT procedures and simple estimates 
of system parameters (e.g., estimates of Eq), especially when capabilities for shaping control 
fields are available. By combining these QCPs, we have developed a form of hybrid quan- 
tum control; estimates of system parameters can be directly incorporated into simulations 
to generate improved quantum operations for information processing and memory. 

Figures 11 and 12 compare the vector-constraint norm ||i7^||2 and objective functional J 
as a function of the optimization iteration for both operations. Overall, ||^/'^||2 increases as 
J decreases, which is consistent with the notion of minimizing ||rf''||2 and J as potentially 
competing control objectives. Even though components of Vj7 that are parallel to all 
gradients Vrji are removed at each iteration, \\'f]^\\2 increases during the optimization for (at 
least) two reasons: (a) eq. (21) removes components of V?7j (where elements rji are nonlinear 
functions of the control) using an iterative linear projection method and (b) convergence of 
the DPC+OCT routine does not depend on ||57'^||2- 

To aid in the comparative analysis of results from OCT, DPC, and the DPC+OCT 
procedures, OC gate-distance data from figs. 3 and 4 are also presented in figs. 9 and 10, 
respectively. Although the OC fields Co{sq', t) all outperform the DPC+OCT fields Ch(eo; 
at these OCs do not have the robustness of the DPs or DPC+OCT fields. To emphasize 
this feature, figs. 13 and 14 present Z^/2 and Z.„ gate distances for Co(2;t), Cd(t), and 
Ch(2;t) controls over a unit interval centered at Eq = 2. For both gates, Co(2;t) is very 
sensitive to variations in e, e.g., when e changes from 2 to 2 ± 0.01 (a change corresponding 
to approximately 1.3 x 10~^ T for the GaAs DQD example), the ^7^/2 gate distance increases 
(approximately) from 10^^ to 10~^, while the Z^^ gate distance increases from 10~^ to 10~^, 
approaching the fault-tolerant threshold. However, for Ch(2;t) for both gates, as e varies 
from Eo, the increase in gate distance is much more gradual. Figures 13 and 14 contain 
some useful information to help understand the benefit of the hybrid DPC+OCT protocol. 
By combining DPC, OCT, a DP that satisfies 57'' = for the initial GrA iteration, and an 
estimate of the value of eo, the gate distance is decreased compared to the gate distance of 
the original DP for the entire unit interval centered at Eq. Depending on the uncertainty 
magnitude of e, this benefit could yield a potentially substantial decrease in the required 
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concatenation/encoding resources necessary for QECCs, which depend on gate errors. 

As a final illustrative example, consider the operation applied to the initial state 
|cr^), implemented with the corresponding controls Co(2;t), Cdit), and Ch(2;t), which are 
applied to an ensemble of systems described by the Hamiltonian in eq. (1) and the interval 
1.5 < e < 2.5 (numerically distributed over 20 equal increments of 0.05 scaled units). The 
target state is \a~) = Z^^la^), where (JxWt) — M'^t)^ ensemble of final states for a 

given control is denoted by {iV'i)}- This state-based example clarifies the gate improvement 
obtained from Ch(2;t), compared to Co(2;t) and Cd{t). To quantify the fidelity of the 
controls, we use the Uhlmann state fidelity for pure states [79, 80]: 

^u(|^i),|^2)):=|(^i|^2)|, (25) 

where iV'i) and \1IJ2) are normalized vectors in "H. For a given control, we denote the resulting 
minimum, maximum, and average state fidelity, and the standard deviation of the fidelity 
of the ensemble as min|^^) J^^, max|^^) J^^, and ctjt^, respectively, which are presented in 
table V for Co(2;t), Cd(t), and Ch(2;)f:). Comparing the respective quantities, Ch(2;t) has 
the largest average and minimum fidelity and the smallest standard deviation of fidelity of 
the ensemble (by nearly a factor of 10). The final-time ensembles are also illustrated in 
fig. 15, which contains a plot of resulting final states for each control, along with the target 
state \<J~), all in the Bloch vector coordinates y and z. Because — 1 < x < —0.995 for all 
final states, it is not included in this figure. Unlike the Bloch vector components corre- 
sponding to the final states produced from Co(2;t) and Cd(t), the Bloch vector components 
produced from Ch(2;t) are tightly distributed around the target state, with most of the 
error distributed uniformly along the z-axis, centered at the target state |cr~). Very similar 
results are obtained for the operation implemented with Co(2;t), Cd(t), and Ch(2;t), 
applied to an ensemble of systems where 1.5 < e < 2.5. 

VIII. CONCLUSIONS AND FUTURE DIRECTIONS 

Combining OCT with the DPC established by Pasini et al. introduces improvements to 
the control of quantum systems for information processing. Given a reasonable characteri- 
zation of the angular frequency e/h of a persistent, but somewhat uncertain rotation about 
the X-axis, a near-optimal fidelity can be achieved for e > 0. Furthermore, the resulting 
DPC+OCT controls exhibit improved robustness to uncertainty in e, compared to the origi- 
nal DPs. The systematic integration of general DPC and control-field shaping methods from 
OCT, therefore, promises considerable improvement over DPC or OCT strategies alone. We 
have provided a quantitative illustration for a logical qubit based on a DQD system, with 
continuous controls that possess reasonable magnitudes [64], based on the scaled-to-SI unit 
mapping. 
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Target operation: 








1 


2 


3 


4 


5 


max Ch 


29.5 


29.5 


29.3 


29.2 


29.4 


29.6 


^(tf;Ch) 


7r/2 


7r/2 


1.5704 


1.5690 


1.5695 


1.5792 




335.5 


334.6 


332.8 


337.4 


356.5 


370.5 


A(Z,/2,C/tf) 


3.65 X 10"*^ 


1.21 X 10-6 


8.23 X 10-6 


1.86 X 10-5 


4.08 X 10-6 


2.55 X 10-6 


ll^'^(if;Ch)||2 


3.21 X 10"*^ 


4.74 X 10"^ 


2.13 X 10-3 


4.21 X lO--'^ 


4.12 X 10-3 


2.45 X 10-3 


7e,/2[Ch, £0,0.5] 


2.84 X 10-6 


9.47 X 10-5 


3.55 X 10-"^ 


6.46 X 10-^ 


5.72 X 10-^ 


9.20 X 10-"^ 



Target operation: Z-,^ 








1 


2 


3 


4 


5 


max Ch 


28.8 


28.8 


28.6 


28.4 


28.4 


28.9 


^(tf;Ch) 


vr 


vr 


3.1411 


3.1392 


3.1343 


3.1320 




264.8 


264.0 


261.8 


259.1 


260.0 


280.7 




2.36 X 10-*^ 


6.57 X lO-'^ 


1.67 X 10-5 


1.74 X 10-5 


2.79 X 10-*^ 


1.10 X 10-6 


ll^'^(if;Ch)||2 


3.55 X 10-9 


9.93 X 10-5 


4.90 X 10-"^ 


1.10 X 10-3 


6.49 X 10-3 


2.10 X 10-2 


7^<^[Ch, £0,0.5] 


1.83 X 10-5 


3.11 X 10-^ 


1.18 X 10-3 


2.59 X 10-3 


4.23 X 10-3 


4.92 X 10-3 



TABLE IV: Performance of the DPC+OCT controls Ch(eo; t) for one-qubit operations. 
Here, max|Ch|, Q-, ^iC*] := ^^^C'^{eQ\t)^t^ A, ||rf''||2, and Jl,^ are the maximum control-field 

amplitude, angle of controlled 2;-axis rotation, control-field fiuence, gate distance, 
constraint vector norm, and gate robustness, respectively, in the corresponding scaled units 

described in section II C. 



We are currently investigating the benefits of these DPC-I-OCT 7r/2- and vr-pulses for 
memory and information processing in the presence of a decohering spin bath. It will 
be useful to determine how these pulses extend spin-echoes and improve general DD and 
dynamically-corrected gate pulse sequences, such as those described in refs. [3, 4, 17, 81- 
85] . Future work involves an exploration of unitary control sensitivity to fluctuations in the 
control field itself (e.g., control noise). Post facto analysis of both the OCT and DPC+OCT 
results presented in this article for the general qubit model suggests that these fluctuations 
may contribute just as significantly to gate errors as corresponding system and environment 
fluctuations. However, our optimization criteria does not include robustness to control-field 
noise; such robustness may be sacrificed in favor of the actual criteria. Given the ubiquity 
of noise in classical controls and quantum-mechanical systems, constructing controls and 
systems that are robust to their own noise is crucial for practical fault-tolerant QC. 

Extensions of the original analysis by Pasini et al. are also being considered. We are 
interested in generalizing their results to include (a) arbitrary angle rotations and axes, (b) 
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Target state: |cr^ ) = Z.„\a^) 













Co(2;t) 


0.996197 


1.0 


0.998871 


1.082 X lO"-'^ 


Cd(i) 


0.999678 


0.999985 


0.999884 


9.398 X 10"^ 


Ch(2;t) 


0.999958 


1.0 


0.999991 


1.129 X 10-^ 



TABLE V: Properties of state fidelity for the transition \a~) = Zj,\a^) driven by controls 
Co(2;t), Cd(t), and Ch{2;t). Here, min|^-) J^u, max|^^) and ajr^ denote the 

minimum, maximum, and average fidelity, and the standard deviation of the fidelity, 
respectively, of the ensemble of final states compared to the target state \cr~). 

closed-system perturbative expansions about any value of e, rather than only e = 0, and (c) 
e as a stochastic time- dependent variable/operator, which is relevant to previous research 
on decoherence control, e.g., [86]. In addition, direct minimization of TZ^ [eq. (13)] or 

) de (26) 

HS/ 

over J-\0), where Cl, C2 G M weight the relative significance of the two norms, are purely 
OCT means to improve robustness to variations in e about any fixed interval [£1,6:2], which 
we are also investigating. 
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Time t (scaled units) 



FIG. 1: (Color online) OC fields producing a Zt^i2 operation for several specific values of Eq 
and tf = 1 scaled unit. All OCs have distances A{Zt^/2, Ut^) < 10~^. The inset presents the 
difference between Co(l;t) and Co(2;t). Although distinct, note that Co(l;t), Co(2;t), and 
Co (3; t) appear nearly indistinguishable in this figure. 
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FIG. 2: (Color online) OC fields producing a operation for several specific values of eo 
and tf = 1 scaled unit. All OCs have distances A(Z.n-, f/tj < 10~^. The inset presents the 

difference between Co(l;t) and Co(2;t). 
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FIG. 3: (Color online) Distance of OCs optimized with particular values of epsilon, 
denoted as Co{eo;t), for the operation [A(Z.n-/2, f/tj < 10~^ for all controls], 
subsequently applied over the interval < £ < 6 (with a resolution of 0.01 scaled units). 
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FIG. 4: (Color online) Distance of OCs optimized with particular values of epsilon, 
denoted as Co{eo]t), for the operation [A(Z^,f/tj) < 10~® for all controls], subsequently 
applied over the interval < e < 6 (with a resolution of 0.01 scaled units). 




FIG. 5: (Color online) Control fields satisfying the DPC in eq. (17), denoted as Cd(t), for 

operations {f]^\\2 < 10"^), from Pasini et al. [15]. 
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FIG. 6: (Color online) Distance of the operations as a function of e (with a resolution 
of 0.01 scaled units) and G {7r/2,7r} for the Landau-Zener model of eq. (1) and control 
fields Cd(t) presented in fig. 5. The inset displays the gate distance for < £ < 1 in greater 

detail, on a linear scale. 
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FIG. 7: (Color online) DPC+OCT fields Ch(£o;^) producing a Zt^i^ operation, optimized 

using Cd(t) as the initial control for all estimates/values of Eq and tf = 1 scaled unit. 
Although distinct, note that Ch(0;t), Ch(l;t), and Ch(2;t) appear nearly indistinguishable 

in this figure. 
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FIG. 8: (Color online) DPC+OCT fields Ch{£o',t) producing a Z.,, operation, optimized 
using Cd(t) as the initial control for all estimates/values of Eq and tf = 1 scaled unit. 
Although distinct, note that Ch(0;t), Ch(l;t), and Ch(2;t) appear nearly indistinguishable 

in this figure. 
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FIG. 9: (Color online) Distance for the Z-j^ji operation as a function of e (with a resolution 
of 0.01 scaled units) for control fields satisfying (a) i/'' = [red dashed line, which is very 
similar to that for Ch(0; t)], (b) ff and A 0, optimized with a specified value of Eq, 
and (c) results from Co{so;t) in section IV (black dashed lines, from fig. 3). 



34 




1 2 3 4 5 6 

e (scaled units) 



FIG. 10: (Color online) Distance for the operation as a function of e (with a resolution 
of 0.01 scaled units) for control fields satisfying (a) i/'' = [red dashed line, which is very 
similar to that for Ch(0;t))], (b) f/"" f« and A 0, optimized with a specified value of Sq, 
and (c) results from Co{sQ]t) in section IV (black dashed lines, from fig. 4). 




Iteration Iteration 

FIG. 11: (Color online) ||57'^||2 and objective functional J for the ^7^/2 operation and each 
value of So considered, as a function of the number of DPC+OCT algorithm iterations. 
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Iteration Iteration 

FIG. 12: (Color online) \\'f]'^\\2 and objective functional JT" for the Z^^ operation and each 
value of considered, as a function of the number of DPC+OCT algorithm iterations. 




e (scaled units) 

FIG. 13: (Color online) Distance for the ^7^/2 operation as a function of e for Co(2; t), 
Cd(t), and Ch(2; t) controls in a unit interval centered at eo = 2 (with a resolution of 0.01 
scaled units). Note that for e > 0.7, A[Z^, U{tf, Cd)] > A[Z^, U{tf, Ch)], i.e., Ch(l; t) 

outperforms Cd(t). 
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FIG. 14: (Color online) Distance for the operation as a function of e for Co (2; t), Cd(t), 
and Ch(2; t) controls in a unit interval centered at eg = 2 (with a resolution of 0.01 scaled 
units). Note that for e > 1.5, A[Z^,U{t{;Cd)] > A[Z„, f/(tf; C^)], i.e., Ch(l;t) outperforms 

Cd(t). 
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FIG. 15: (Color online) Final state in the Bloch vector coordinates y and z for the 
operation, implemented with the corresponding controls Co(2;t), C^it), and Ch(2;t) 
applied to an ensemble of systems described by the Hamiltonian in eq. (1) and the interval 
1.5 < e < 2.5 (distributed over 20 equal increments of 0.05 scaled units). The target state 
is W^) = Z^|cr+). Because —1 < x < —0.995 for all final states, it is not included in this 

figure. 
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